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It is commonly believed that the quantum wave-function of a very massive particle with small 
initial uncertainties in both position and velocity (consistent with the uncertainty relation) will 
spread very slowly, so that the limiting dynamics is classical Newtonian and, in particular, deter¬ 
ministic. This argument assumes, however, that the classical motions for given initial data must be 
unique. It was discovered about two decades ago in the field of fluid turbulence that non-uniqueness 
due to “roughness” of the (non-random) advecting velocity field may lead to stochastic motion of 
classical particles. Vanishingly small random perturbations are explosively magnified by turbulent 
Richardson diffusion in a “nearly rough” velocity field, so that motion remains stochastic as the 
noise disappears. This classical spontaneous stochasticity effect has great importance for magnetic 
reconnection and dynamo in astrophysics, for geophysical mixing problems, and for everyday mani¬ 
festations of fluid turbulence. Simple analogies between stochastic particle motion in turbulence and 
quantum evolution (Wiener vs. Feynman path-integrals, Fokker-Planck vs. Schrodinger equations) 
suggest that there should also be quantum spontaneous stochasticity (QSS) driven by quantum fluc¬ 
tuations. We verify this idea for ID models of a quantum particle in a repulsive potential which is 
“nearly rough”, with V{x) ~ —C\x\^'^°‘ at distances \x\ ^ £ for some physical UV cut-off I, and for 
an initial Gaussian minimum-uncertainty wave-packet centered at 0. We first consider the WKB 
limit with h/m —>■ 0, then position-spread u —>■ 0. As in critical lattice systems, we show that the 
“infrared” limit of time t —>■ cxa is equivalent to the “continuum” limit l,a ^ 0, with p, — Ifa fixed. 
The corresponding scaling limit of the position density is non-deterministic, with equal probabili¬ 
ties of the particle following either of two, non-unique classical solutions and with super-ballistic, 
Richardson-like spreading of the wave-packet. The splitting of the wave-packet into two peaks oc¬ 
curs in a very short time t ~ rc(/i)(f!^^“/C)^'^^ for small 1. We derive the asymptotics of for 
small and large p. We also consider the scattering off the rough potential of minimum-uncertainty 
wave-packets with non-vanishing expectations {x), {p) and find that QSS occurs here too for careful 
fine-tuning of {x), (p), a, 1. Beyond WKB, we also consider other semi-classical limits by a nu¬ 
merical solution of the Schrodinger equation, where we observe QSS both in position-space and in 
momentum-space. Although the wave-function remains split into two widely separated branches in 
the classical limit, rapid phase oscillations within each branch prevent any coherent superposition. 
Thus, indeterminism persists in the classical limit for rough potentials, but not quantum entan¬ 
glement. These are unambiguous predictions of quantum mechanics which should be observable in 
laboratory experiments. We consider possible realizations by AMO, ultra-cold neutrons, and optical 
analogue experiments. We also briefly contrast QSS with Anderson localization. 


I. INTRODUCTION 

The determinism of Newtonian particle dynamics is 
generally taken for granted. However, the uniqueness of 
solutions of classical initial-value problems for ordinary 
differential equations (ODE’s) requires mathematical as¬ 
sumptions on smoothness of the velocity vector field, and 
there are simple counterexamples in every textbook on 
ODE’s when these assumptions are not satisfied (e.g. see 

[1] , section II.5.) Such examples have engendered philo¬ 
sophical ruminations on acausality in classical dynamics 

[2] . This subject entered the realm of physics in the fun¬ 
damental work of Bernard, Gaw§dzki and Kupiainen on 
particle advection in turbulent fluid flow [3], which pre¬ 
dicted that Lagrangian particle trajectories must become 
non-unique and indeterministic for an exactly prescribed 
fluid velocity field and an exactly prescribed initial par¬ 
ticle location, as a consequence of the “roughness” asso¬ 
ciated to a Kolmogorov energy spectrum. Those authors 
pointed out that stochasticity in solutions of determinis¬ 


tic ODE’s should occur in the situation physically real¬ 
ized by fluid turbulence with “roughness” of the velocity 
field only over a finite inertial-range of scales—and thus 
smooth velocity fields at sufficiently small scales—when 
coupled with vanishingly small random perturbations to 
fluid trajectories. This phenomenon was later dubbed 
spontaneous stochasticity [4]. It is, in retrospect, a di¬ 
rect consequence of the pioneering 1926 work of L. F. 
Richardson on turbulent particle dispersion [5], but long 
unappreciated. This phenomenon has profound conse¬ 
quences in astrophysics, for example, where spontaneous 
stochasticity of plasma fluid motions breaks the ideal con¬ 
straint of magnetic flux-freezing even in the limit of infi¬ 
nite conductivity and allows magnetic reconnection to be 
fast [6-9]. Analogous breaking of the Kelvin-Helmholtz 
theorems on vortex-line motions occurs in hydrodynamic 
turbulence [10], with consequences for turbulent mixing 
of vorticity in geophysical and engineering flows. 

The empirical verification of spontaneous stochastic¬ 
ity effects in fluid and plasma turbulence currently rests 
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upon numerical simulations of the Navier-Stokes and 
magnetohydrodynamics equations [9, 11, 12]. The phe¬ 
nomenon has not yet been seen in controlled laboratory 
experiments of turbulent flows [13], presumably because 
of the technical difficulties in tracking advected parti¬ 
cle pairs that start sufficiently close together. The best 
prospect for fluid turbulence experiments is to make a 
bigger flow apparatus with a larger Reynolds number. In 
this way, particles even at relatively large initial sepa¬ 
rations have sufhcient time within the inertial-range to 
exhibit the effects. An alternative idea that we explore 
here is instead to go smaller and to search for analogous 
effects in the quantum world. There is in fact intrinsic 
interest in the question of how a quantum system will be¬ 
have whose corresponding classical dynamics is (nearly) 
non-unique. As we shall see, this situation leads to quan¬ 
tum randomness surviving in the classical limit where one 
normally expects determinism to prevail. Before demon¬ 
strating this, we first review briefly the phenomenon of 
classical spontaneous stochasticity. 

II. BACKGROUND 

The results of Richardson’s theory of particle disper¬ 
sion in turbulent flow can be most easily understood from 
a toy model for the growth of separation r between a pair 
of particles advected by the flow: 

dr/dt = C'(er)l/^ (1) 

where the righthand side represents the relative velocity 
6u{r) of the particles at separation r. The scaling with 
exponent 1/3 is consistent with the statistical theory of 
Kolmogorov [14]. Solving by separation of variables gives 
the exact solution 

which exhibits the Richardson scaling r^(t) cic et^ [5] for 
t 3> . More surprisingly, however, 

/9 

r{t) ^ > 0 

as ro —> 0. Thus, two particles at identically the same 
initial point, with tq = 0, end up a finite distance apart 
at any positive time! One’s natural expectation would 
be instead that r{t) = 0 and that the two particles would 
move together for all time. In fact, these are just two 
examples of infinitely-many solutions of the initial value 
problem for equation (1) with rg = 0. The most general 
solution has, for any waiting time r G [0, oo], r{t) = 0 for 

t < T and then r{t) = for t > t. 

The above toy problem is the standard example of gen¬ 
eral theorems on the non-uniqueness of solutions of the 
initial-value problem for an ODE: 

dx/dt = u{x,t), x(to) = xg. (2) 


If the advecting velocity field satisfies a spatial Holder 
continuity condition 

lu(x-I-r) - u(x)l < Clr]" (3) 

for a = 1 (a so-called Lipschitz condition) then solu¬ 
tions are unique, but examples like the previous one with 
a = 1/3 show that solutions need not be unique for 
0 < a < 1 [1]. The latter situation is the expected one for 
turbulence. At the level of individual flow realizations, it 
was shown by Onsager [15, 16] that the observed energy 
dissipation anomaly for turbulent flow requires that in 
the limit of infinite Reynolds number the fluid velocity 
field must have a Holder exponent a < 1/3 and, in par¬ 
ticular, u cannot remain space-differentiable as Re —)■ oo. 

Real fluid turbulence has, of course, a finite Reynolds 
number albeit possibly very large. Thus, the velocity 
field is ultimately smooth and satisfies a bound of the 
form (3) with a = 1 for jrj < ry, a small length-scale set 
by viscosity v. For example, in the Kolmogorov theory, 
r] = Thus, particle trajectories are unique. 

However, it was realized in [3] that the non-uniqueness 
resurfaces for the fluid trajectories with small random 
perturbations to the dynamics or to the initial data. That 
paper considered the specific problem of stochastic par¬ 
ticle trajectories satisfying a Langevin equation 

dx/dt = u^(x,t)-I-j7(t), x(tg) = Xg-I-p I (4) 

with u^(x, t) the turbulent fluid velocity field, with f){t) 
a vector white-noise, and with ^ any random vector with 
zero mean and finite variance. The subscript ly stands 
for the (kinematic) viscosity of the fluid, which is a finite, 
positive value. This stochastic evolution is relevant to the 
problem of an advected passive scalar 9 (e.g. temperature 
or dye concentration) satisfying 

dt9{x, t) Uy (x, t)-Vd(x, t) = kA9{x, t). (5) 

In fact, if p^(x', t'jx, t) is the transition probability den¬ 
sity for the stochastic evolution, then 

9{x,t) = J d^xo 9{xo,to) Pi.(xg,tglx,t), to < t. 

It is convenient to express the transition probability by 
a Wiener path-integral [3, 11, 17] 

Pi.(x,tlxo,tg) = /" d‘’*(x-|(t)) 

•^^(to)=xo 

xexp(-^^ dsli(s) - Uj,(|(s),s)l2) (6) 

If the velocity field u remains smooth (i.e. if the viscosity 
n is held fixed), then a steepest descent argument on the 
path-integral shows that as k —>■ 0 

p^(x,tlxg,tg) ^ d‘’*(x- Xt(xg)) 

where xt(xg) is the unique solution of the initial-value 
problem (2). However, if the turbulent velocity field be¬ 
comes only Holder continuous, as when taking the limit 
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K —)■ 0 together, then it was shown in [3] that a non¬ 
trivial limit p*(x, t|xo,to) may result. Thus, the particle 
trajectories remain random as the external noise is re¬ 
moved and the deterministic problem (2) is recovered. 

There is a close analogy of this “spontaneous stochas- 
ticity” with the zero-temperature phase transition in 
the ID Ising model [18]. We may write the ID Ising 
model in finite-volume [—fV, ...,fV] as a Gibbs distribu¬ 
tion PmIo] = ^ exp for Hamiltonian 

i = -N 

<^-N ^ +1 

with boundary condition CT-at = -1-1 and free bound¬ 
ary condition on ctat. To emphasize the analogy with 
the path-integral (6), we have rewritten the Ising model 
Hamiltonian as a sum of squares and imposed a bound¬ 
ary condition only at the left endpoint. In the zero- 
temperature limit at fixed iV, 

N 

Pn[o'] —> 3'S T —> 0, 

i=-N 

the unique ground-state with a-N = +1. On the other 
hand, if the infinite-volume limit TV —>■ oo is taken first 
before taking the zero-temperature limit, then 

T’oo[ct]2 + 2 T ^ 0, 

i i 

the symmetric mixture of inhnite-volume ground-states. 
This is the well-known zero-temperature phase transition 
in the ID Ising model. The phenomenon discovered in 
[3] for non-smooth dynamical systems is exactly analo¬ 
gous, with the viscosity v corresponding to IV, the noise 
K to temperature T, and the non-unique solutions of the 
initial-value problem (2) corresponding to the non-unique 
ground states of the infinite-volume Ising model. 

The results of [3] were obtained for a soluble turbu¬ 
lence model, the Kraichnan model [19, 20], in which the 
fluid velocity held is replaced by a “synthetic turbulence” 
consisting of a Gaussian random held which is delta- 
correlated in time. The results on spontaneous stochas- 
ticity hold with probability one for every realization of 
the Gaussian velocity held. These results in the Kraich¬ 
nan model have been conhrmed and further elaborated 
in many subsequent works [4, 21-25]. The current evi¬ 
dence of these effects in huid turbulence described by the 
Navier-Stokes equations consists of numerical studies of 
mean-square pair-dispersion (r^(t)), where the average is 
over the initial space position x of the pair in a single how 
realization. These studies conhrm the Richardson predic¬ 
tion (r^(t)) oc and, even more importantly, provide 
evidence that the pair-dispersion for particles evolving 
under the stochastic evolution (4) becomes independent 
of K after a short time of order [11], or, simi¬ 

larly, that the pair-dispersion for particles evolving un¬ 
der the deterministic equation (1) becomes independent 


of the initial separation rp after a time of order 
[12]. Indirect evidence for spontaneous stochasticity is 
also afforded by experimental observation of anomalous 
dissipation of passive scalars in huid turbulent hows [26]. 

There is a close similarity of the path-integral (6) 
for the transition probability of the stochastic advection 
problem and the Feynman path-integral for the tran¬ 
sition amplitude in quantum theory [27]. Consider a 
non-relativistic quantum-mechanical particle of mass m 
and electric charge q moving in an electric held E = 
— Vd" — (l/c)9tA and magnetic held B = VxA gov¬ 
erned by the Schrodinger equation 

ihdtip = — (—ih'V — 

2m V c / 

Then Feynman represented the transition amplitudes by 
the formula 


(x,tlxo,0) = 




<(0)=Xo 


Pxexp ( — / (is L(x( 5), x(s), s) 

^ .In 


where the classical Lagrangian is 

T(x, X, t) = -mlxj^ -I- - A(x, t) • X — q^{x, t). 

2 c 

Not only is Feynman’s formula analogous to (6), but it 
can be used to derive (6) ([11], Appendix). This close 
relation naturally raises the question whether quantum- 
spontaneous stochasticity effects can arise in the classical 
limit, given formally by 7i —>■ 0. When E, B are Lipschitz, 
then the classical equations of motion 


mx = q 


E(x,t) 


1 

c 


X X B(x, t) 


have unique solutions and the stationary phase argument 
of Feynman [27] yields classical dynamics for h —>■ 0. How¬ 
ever, if the electromagnetic helds are non-Lipschitz, then 
quantum randomness could possibly persist in the clas¬ 
sical limit. 

Similar questions have been raised previously in the 
literature. An interesting recent mathematical study 
[28, 29] has considered the classical limit for the sim¬ 
plest ID Hamiltonian system whose classical solutions 
are non-unique. See eqs. (7)-(9) below. In this example a 
power-law potential leads to a force field which is Holder 
continuous with an exponent a < 1. It was proved in [28] 
that the Wigner function for the quantum particle con¬ 
verges in the formal classical limit —)■ 0 to a non-trivial 
probability measure, so that indeterminacy persists in 
the limit. This mathematical result assumes, however, 
an exact power-law potential down to zero-length scale, 
i.e. to the Planck length and below! The physical va¬ 
lidity of the nonrelativistic Schrodinger equation with a 
classical potential will break down well before this scale, 
e.g. quantum electrodynamic effects such as particle pair- 
production will become relevant at the Compton wave¬ 
length h/mc. In any case, an exact power-law down to 
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zero length-scales should not be necessary to obtain such 
effects, just as [3] showed for the problem of turbulent 
advection. The main aim of our work will be to in¬ 
troduce physical short-distance cut-offs into the model 
of [28] and to show that quantum spontaneous stochas- 
ticity effects exist in experimentally accessible regimes. 
The need for such a short-distance cut-off was realized 
in previous work of Bliimel [30], who studied a related 
phenomenon of “quantum ray-splitting” in discontinu¬ 
ous step-potentials (see also [31, 32]). In this case also, 
quantum randomness persists in the classical limit of van¬ 
ishing de Broglie wavelength A as the short-distance cut¬ 
off £ is taken simultaneously to zero. Such step-potentials 
lead, however, to very singular delta-function forces. The 
innovation of our work will be to show that such effects 
can occur with much smoother force fields which, even as 
the cut-off i is removed, remain bounded and continuous. 


III. MODELS 


which reduces the system of two equations to a single 
first-order ODE for |x(t)|: 


d 

dt 


|x| = ± 



\x\^+-+Ho, 


x{to) = Xo 


( 10 ) 


with ± representing sign(a: • v) and Hq = ^Vq + V{xo)- 
For any Hq < 0, one can write Hq = — where 

xq satisfying [xo] < |a;o| is the position of a particle in¬ 
stantaneously at rest with the same energy. Separating 
variables yields [33] an implicit solution in terms of hy¬ 
pergeometric functions ( [34], Ch.l5; [35], Ch.II): 




2’ 2{l + a)’ 2{l + a)’ 


1+a 


2’ 2(l + a)’2(l + a)’ 


l+aN 




2C 
\ a 


1/2 


(t-to). 


( 11 ) 


Just as in [28], we consider the model of a single parti¬ 
cle moving in a repulsive, scale-invariant cusp (or conic) 
potential energy field 

E(x) = -^Ni+“, 0<a<l. (7) 

1 -I- a 

In fact, we take V to be the potential energy per mass m 
of the particle, so the constant C > 0 absorbs a factor of 
1 /m. The classical dynamics is thus given by the ODE 

X = V, V = C'|a;|“sign(a;), (8) 

which are Hamiltonian with H = -I- V{x). The flow 

vector field U(a;,u) = (u, C'|a;|“sign(a;)) is infinitely dif¬ 
ferentiable except at points of the form (x,v) = (0,uo), 
where it is only Holder continuous with exponent a. In 
particular. 


|U(r,0)-U(0,0)| = qr|“, 


at point (x,u) = (0,0). It is not hard to show that the 
solutions of the initial-value problem are unique at all 
starting points except those that lead to solutions pass¬ 
ing through (x,v) = (0,0), where there are infinitely 
many solutions. The two extremal solutions starting at 
(xo,vo) = (0,0) for initial time to = 0 (which reach the 
farthest distance at each later time t) are 


x±(t) = ± 


1 

2 


(I-a) 




2/(l-a) 


(9) 


There are infinitely many other solutions for this same 
initial condition in which the particle waits at a; = 0 for 
a time r € [0, oo], before moving off in the same manner. 

The general solution of the initial-value problem for (8) 
can be obtained by quadrature using energy conservation, 


If sign(xo • Vo) < 0, then the — sign holds initially 
and I a; I decreases to [Sq], whereupon the velocity re¬ 
verses direction, |a;| increases and the -|- sign holds. For 
ffo > 0 there is a similar formula obtained by writing 
^0 = but with the -I- sign in the main ar¬ 

gument of the hypergeometric function changed to a — 
sign. In this case, v never changes sign, but |a;| may pass 
through 0 with a corresponding of change of overall sign 
on the righthand side from — to -I-. It is easy to see from 
the exact solution (II) that, for all choices of initial con¬ 
ditions (xo,vo), the ratio x(t)/x±(t) —)■ I as t —)■ oo for 
one choice of ±. This “forgetting” of the initial data is 
the essence of spontaneous stochasticity. 

A physically realizable potential will however exhibit 
a power-law scaling as in (7) only over a limited range of 
x-values, with an “inner” or short-distance cutoff I and 
an “outer” or large-distance cut-off L. We shall not in¬ 
troduce an explicit cutoff L at large-distances, but we 
can use our results with L = oo to estimate the practical 
effects of such a cut-off. We introduce two concrete mod¬ 
els with a short-distance cut-off i. The first is a globally 
analytic potential given by a confluent hypergeometric or 
Kummer function ([34], Ch.I3; [35], Ch.VI; [36], Ch.IV): 


Vt,{x) = - 


^gi+^ 2(i+-)/^r((2 + a)/2) 

(I -I- a)v^ 

x.n(-l(l + a),l;-A 


( 12 ) 


We shall refer to this as the “Kummer potential”. It can 
be obtained by convolving the power-law potential in (7) 
with a Gaussian smoothing kernel exp(—x^/2£^)/-\/27r^, 
and thus becomes indistinguishable from the power-law 
for I a; I £. The second model potential is chosen to 

coincide exactly with the power-law for |a;| > I but is 
only piecewise smooth with continuous first-derivatives 
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at |x| = £, given by the elementary formulas: 



We shall refer to this as the “first-order spliced poten¬ 
tial”. It has less generic behavior because of its lower 
degree of spatial smoothness but it is more tractable for 
analytic calculations. 

The main problem that we shall consider is the quan¬ 
tum evolution in these potentials of an initial minimum- 
uncertainty Gaussian wave-packet 

(■ ** S’*'**) 

with mean {x) and uncertainty a in position-space and 
mean {p) and uncertainty in momentum-space. 

There is a clear analogy with the problem of turbulent 
transport, which we make explicit in the table below: 


Turbulence-Quantum Analogy 


Turbulent Advection 

Quantum Mechanics 

Fokker-Planck equation 
dtP = K,Ap — V-(up) 

Schrodinger equation 

-t 

Velocity integral length L 

Outer length L of potential 

Dissipation length rj 

Inner length i of potential 

Scalar diffusivity k 

Quantum of circulation h/m 

Position uncertainty p 

Wavepacket spread a 

Reynolds number Re = (LjrkA Scale ratio Lji of potential 


stressed that the WKB limit is not necessary for quan¬ 
tum spontaneous stochasticity but it is probably one of 
the easiest regimes to study experimentally and, also, it 
has the conceptual advantage that the relation between 
the quantum and classical dynamics is transparent. 

We apply the standard time-dependent WKB approx¬ 
imation [41-43], which corresponds to making the ansatz 

= exp(iQ{x,t)/%) 

for the solution of the Schrodinger equation, together 
with the asymptotic expansion 

0(x,t) = S(x,t) + (h/i)Si(x,t) + (l%/i)‘^S 2 {x,t) -\ -. 

Dropping quadratic terms and higher, this gives 

'ip{x,t) = a/ p{x, t) exp{iS{x, t)/h) 

with p{x,t) = exp(2S'i(a:,t)), and plugging into the 
Schrodinger equation yields the Hamilton-Jacobi equa¬ 
tion (or eikonal equation) for the classical action S{x,t) : 

dtS{x,t) + + Vix) = 0, 

2m 

together with the conservation (transport) equation for 
the position probability density p{x,t) : 

dtp{x, t) + V ■ {v{x, t)p(x, t)) = 0 

with v{x,t) = VS{x,t)/m the classical velocity. These 
equations are most easily solved by method of character¬ 
istics, using Hamilton’s equations of motion 


Note that the Fokker-Planck equation for the probability 
distribution p(x, t) of an advected particle is the adjoint 
to the passive scalar equation (5). Similar analogies be¬ 
tween high-Reynolds-number turbulence and the classi¬ 
cal limit of Schrodinger evolution have been pointed out 
already some time ago (e.g. [37], section 6), but little 
pursued up until now. 


x=p/m, p=—W{x) (15) 

with the initial conditions a;(0) = xq andp(O) = VS'o(a:o)- 
Calling the resulting solutions xt(xo) and pt{xo), it fol¬ 
lows that S{xt{xo),t) = St{xo) with 

^((xo) = S'o(a;o) -f J ds [^pKxq) - V{xs{xo))]. (16) 


IV. WKB ANALYSIS 


The probability conservation equation can likewise be in¬ 
tegrated along characteristics as 


We shall first analyze our model problem in the limit 
where the de Broglie wavelength A = h/mv of the par¬ 
ticle is much smaller than the inner length £ of the po¬ 
tential. This is the usual semi-classical or WKB limit. 
In the turbulence-quantum analogy sketched above, this 
corresponds to the limit of very large Prandtl number 
Pr = vjK. As discussed elsewhere ([11], section IV) 
classical spontaneous stochasticity effects are expected to 
be observable for Lagrangian particle dispersion in high- 
Prandtl number turbulent flows at sufficiently long times. 
This is an importance astrophysical regime appearing in 
the interstellar [38] and intracluster [39] media and rel¬ 
evant to problems of star formation [40]. It should be 


p{xt(xQ),t) 


Po(a:o) 

Jt{xo) 


where we have introduced the Jacobian 


(17) 


Jtixo) = dxt{xo)/dxo = x^{xt)). 

To calculate the latter it is most convenient to differen¬ 
tiate Hamilton’s equations to obtain a set of ODE’s 

j = AT, k = -V''{x)J (18) 

for Jt{xk and Kt{xo) =k{xo). 
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We apply these standard results to our model prob¬ 
lems. Our scaling of V corresponds to setting m = 1. The 
initial Gaussian wave-packet (14) corresponds to taking 


So{x) = {p)-x, poix) = 


jx - jx))^ 

2ct2 


In that case, we have initial conditions in the ODE’s 

a:(0) = xo, p(0) = (p), J(0) = 1, K{0) = 0. 

There are no caustics for this problem, so that the WKB 
approximation is globally valid in time. As usual, the 
WKB approximation is locally valid in space when 

9 / \ ( dS 1 

a ■I’"'} 

where A(a:, t) = h/p{x, t) is a local de Broglie wavelength 
with p{x,t) = dS{x,t)/dx. E.g. see [41]. This corre¬ 
sponds to our starting assumption X <^ £ when the clas¬ 
sical action S{x, t) develops no smaller length-scales than 
in the potential V{x). Note finally that the Schrddinger 
equation divided by m is 

i — dt'ip = f—^ + Vip (19) 

m 2 \m) 

where V/m —> V now represents in reality the potential 
energy per mass. Thus, holding V fixed, the only de¬ 
pendence of % is through the ratio fi/m and one may, by 
varying m, study the problem formally with m = 1 and 
varying h. (Of course, the dependence of E on m must 
be remembered in comparing results with experiment.) 


A. Scaling Limit 

We shall begin by studying in this and the following 
section the symmetrical problem with {x) = (p) = 0. We 
are interested here primarily in the probability density 
p{x, t) for the location of the particle. Since in the WKB 
limit there is no dependence of p{x, t) upon Planck’s con¬ 
stant ?i, it contains only dimensional parameters cr, t and 
C, in addition to cc, t. It is natural to introduce a dimen¬ 
sionless scaled position 

X = x/x+{t). 

Then simple dimensional analysis gives 

p{x,t]a,£,C) = -^j-)_p(a;/a;+(<),r;p) 

where p is the probability density of x and 

T = p = ila. 

This result tells us that while keeping the ratio p, = Ija 
fixed, the limits t —)■ oo and cr —> 0 are identical. Quan¬ 
tum spontaneous stochasticity will be present if 

p{x, r; p) ^S{x -f 1) -f id(x - 1), (20) 


or, equivalently, if at fixed time t 

p{x,t]a,£,C) —^ ^S{x - x-{t)) + ^5{x - x+{t)) 

in the limit a,t Such a limit was ob¬ 

tained in [28] for the exact power-law potential. We now 
establish it here for the model with UV cut-off £. 

To prove this result, we smear p with a smooth test 
function ip{x), giving 

Using the WKB formula (17), this can furthermore be 
written as 

J dx <p{x)p{x,t) = J dxo p{xt{xo)/x+{t))po{xo), 

where we used the definition Jt = dxt/dxQ to change the 
integration variable to xq. The essence of the result is 
now to show that 

lim xt{xQ)/x+{t) = ±1, (21) 

t—^OO 

with ± = sign(a:o). In fact, this yields a more general 
result than (20), for any choice of po(x) which contains 
no delta function part cx S(x). In that case, we can write 

/ dx (p(x)p(x,t) = / dxo (p(xt(xo)/x+(t))po(xo). 

J J —OO 

-f / dxo ip{xtixo)/x+(t))poixo), 

Jo 

and apply dominated convergence to infer that 

lim / dx ip{x)p{x,T) = ip{—l)p- + ip{l)p-, 

T^OO J 


/ O /•OO 

dxQ Po{xq), p+= dxo po{xo)- 

-OO J 0 

This gives as a special case the result (20) for any sym¬ 
metrical choice of po with p_ = = 1/2. 

We now demonstrate the crucial result (21) for the 
first-order spliced potential. Here we can develop explicit 
solutions for Xt and Vt- Eor xq < ^ it is straightforward 
to show that 


Xt{xo) = Xo cosh 


£(l-a)/2 


when t < ti(xo) with the definition 

£(l-a)l2 

tt(.xo) = ^^ 1 ^ arccosh(£/a:o), 

at which time Xt = £■ Thereafter one can use (11) to 
write 




I - I 1 + CK 

PEO 


) ’ \xt\ 
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FIG. 1 WKB position-space probability densities for ^ = 1 and dimensionless times r = 0.0, 2.5, 5.0, 7.5. Solid is for the 
Kummer potential, dashed for the spliced potential, both with a = 1/3. The two asterisks on the abscissa indicate the 
locations of the extremal classical solutions at the particular time of each snapshot. For movie, see [33]. 


_^(i-a)/2 ^ /I _ (!-«) . 

^2’ 2(1 -ha)’ 


1 -|- 3q: 






= 2(l-«) 


2(1-ha) 
nn \ 1/2 

^ (t-Mxo)). 


(23) 


for t > tiixo) where it follows directly from the definition 
of the spliced potential (13) that 


|Sor+“ = Q(1 - a) + i(l + a) , 


with xq as an “effective initial position.” On the other 
hand, if xq > I, then the particle starts in the region 
where the cusp potential is exact and for all times t > 0 
the result (11) gives 


Q, 


(1 — a) 1 -h 3a 
2(l-ha)’ 2(l-ha)’"^ 


Xt 


(1 — a) 1 -h 3a 
2(1-ha)’ 2(1-ha) 
2(7 


t. 


1 + a, 




(24) 


Here we have used the fact that xq = xq in this region. 
To prove (21) we use these exact formulas and the fact 
that, for (p) = 0, solutions Xt(xo) preserve the sign of 
Xq. Since the righthand sides of both (23) and (24) are 


asymptotic to (x+(t))^^““^/^ as t —)■ oo, one gets easily 
that Xt{xo)/x+(t) sign(xo). This is also true for the 
fixed point xq = 0, with the convention that sign(O) = 0. 

The above results should carry over qualitatively to ev¬ 
ery potential which has the power-law form (7) for some 
long range of x, £ <S1 |x| ^ L, but smoother behavior at 
X i preserving the repulsive character of the poten¬ 
tial. For Xq near 0 there will be an initial period where 
the classical particle position xt(xo) moves away from 
the origin exponentially quickly, until the particle enters 
the region of power-law behavior of the potential. The 
particle will then move outward “explosively”, approach¬ 
ing one of the extremal solutions x±(t) at long times for 
every initial condition xq- We illustrate this with WKB 
results on the position probability density (PDF) for both 
Kummer and spliced potentials. The default values will 
be a = 1/3 and p = 1 in these and all later figures, 
unless stated otherwise. The results were obtained by 
numerically integrating the ODE’s (15),(18) for logarith¬ 
mically evenly spaced initial xq = 10*/'^^ and integers 
i = NsNf,...,NbNp with Nf,Ns < Nb chosen to pro¬ 
vide suitable resolution and then using formula (17) to 
evaluate the position PDF. Because of the symmetry of 
the problem, we could consider only xq > 0. The con¬ 
vergence of the results was tested by varying Nf,Ns,Nb 
and also by verifying that the total probability was unity 
to within at least a percent (and usually much closer). 


























</a:+(t) = 4.4 X 



llx+(t) = 4.35 X IQ-' 



Ijxxit) = 4.35 X 10-^ 



f/a;+(t) = 4.35 x lO^' 



- 0.5 0 0.5 

a:/a:+(t) 


(d) 


FIG. 2 WKB position-space probability densities at a fixed time t for /x = 1 and = 4.35 x 10 10 10 10 

Notations as in Fig. 1. For movie, see [33]. 


For details of the numerical methods, see [33]. Fig. 1 
shows the time evolution of the position PDF’s and their 
splitting from an initial unimodal function into bimodal 
densities at later times. The peaks follow rather closely 
the two extremal classical solutions, indicated by aster¬ 
isks on the abscissa. It is notable that there are jump- 
discontinuities in the PDF’s for the spliced potential, 
which arise from similar jumps in the Jacobians Jt{xo) for 
that potential, due to its lower smoothness. See [33] for 
details. However, qualitatively the PDF evolutions for 
the first-order spliced potential and the smoother Rum¬ 
mer potential are quite similar. In Fig. 2 we illustrate 
instead the convergence of the PDF’s at fixed time t to 
a pair of delta functions for £ = cr <C x+{t), which is the 
key signature of QSS. Again the behaviors for the spliced 
potential and the Rummer potential are quite similar. 

These results have interesting implications for the 
position-space spreading of the initial wavepacket, de¬ 
fined as usual by 

(Aa:(t))^ = J dx (x — (x))^|^(a;, t)p. (25) 


Using our WRB results, we find that asymptotically at 


long times t “/C) 


1/2 


(Ax(t))^ 




the spreading is given by 
1/2 l4/(l-“) 


2C 
I -I- a 


(26) 


This is “Richardson-like” not only in that it is a super- 
ballistic power-law but, more importantly, in that it is 
independent of both h and cr, and thus non-vanishing 
in the mathematical limit of h, a going to zero. This 
should be contrasted with the standard result for ballistic 
position-space spreading of the same initial wavepacket 
under free-particle dynamics: 

(Ax(t))^ = a'^ + {ht/2ma)‘^. 

In this case, the spreading vanishes in the limit h/m —>■ 0 
first (WRB limit) and then cr —0. In fact, this is the 
usual textbook explanation for classical behavior of very 
massive particles. It should be remembered that the con¬ 
stant C in (26) is inversely proportional to m, so that the 
spreading in the cusp potential is also slower for large m, 
but one can always observe a large spread by waiting a 
time proportional to m}/'^ but independent of Ti. 

It is also instructive to compare with the spreading of 
the Gaussian wave-packet in the power-law potential for 
a = 1, that is, the inverted oscillator potential, 

U(x) = -icx2, 

which has been very well-studied [44, 45]. The spreading 
here is: 

(Ax(t))2 = cr^ cosh2(C'i/2t)-b (^ 2mC^/^a ) 
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This result is easily obtained by analytic continuation 
of the well-known result for the harmonic oscillator to 
imaginary frequency. If one takes first the WKB limit 
/i/m — >■ 0 (or h/y/m —>■ 0, since C (x 1/m) and then 
cr —> 0, the spreading vanishes just as for the free par¬ 
ticle. Furthermore, the wavepacket remains unimodal 
and even Gaussian in the WKB limit since the Jacobian 
Jtixo) = cosh(C^/^<) is independent of xq. In fact, this 
is exactly true for the full Schrodinger evolution as well 
[44], and no maxima of at non-zero x ever de¬ 

velop in that case. This example shows that the novel 
features associated to the cusp potential are due to its 
“roughness” or non-smoothness and not merely due to 
strong repulsion from the origin. 

We have ignored so far the effect of a large-scale cut¬ 
off L to the power-law range of the potential. However, 
at a time of the order of the maxima of 

p = ji/P at x±{t) will reach \x\ ~ L and the effects 
of that cut-off will begin to be felt. At that time the 
power-law spreading (26) will terminate. The behavior 
subsequently will depend upon the specific details of the 
potential at length scales |a;| > L. 

B. Splitting Time 

The results of the previous section demonstrate that 
indeterminacy survives as ct, f' —)■ 0 or as t —)■ oo in the 
WKB limit. However, in order to design possible experi¬ 
ments, it is necessary to know exactly how long one must 
wait to observe such effects. The simplest measure of this 
is the time of first appearance of new local maxima of p at 
non-zero values of x, i.e. when the wave-function “splits” 
into multiple peaks. By our previous dimensional anal¬ 
ysis, this will occur at a dimensionless time Tc{p) which 
will depend upon the parameter p, which we refer to as 
the splitting time. We can think of the position x of the 
new maximum as the “order parameter” for transition to 
QSS at the critical time Tc{p). All local extrema of the 
WKB formula p = po/J will occur at x = Xt{xo) for xq 
satisfying (po/J)' = 0, at points where p is differentiable. 
This gives the condition for a local extremum as 

J((xo)/Jt(xo) = Po(xo)/po(a;o), 

or as the vanishing of the Wronskian IF(Ji,po)- 

Because of the symmetry of the classical dynamics 
around the origin J((0) = 0 and because of the symme¬ 
try of the initial wavepacket Po(0) = 0 also, so that there 
is an extremum of p{x, t) at x = 0 for all times. For 
the generic case of a symmetric, smooth potential the 
appearance of the new maxima should appear by a (su¬ 
percritical) pitchfork bifurcation at x = 0, with the ori¬ 
gin changing from maximum to minimum and two new 
maxima appearing to either side. We shall verify this 
expectation for the Kummer potential. Of course, the 
pitchfork bifurcation is unstable to any slight asymme¬ 
try, in either the potential or initial wavepacket, and will 
usually be modified then to a saddle-node bifurcation in 


which a new maximum-minimum pair appear near (but 
not exactly at) x = 0. The observable consequences are 
hardly distinguishable from the pitchfork bifurcation for 
slight asymmetry, so that we analyze only the exactly 
symmetrical problem here. The first-order spliced po¬ 
tential, although symmetric, exhibits also a saddle-node 
bifurcation [33]. In that case the Jacobian Jt(xo) is con¬ 
stant in xo very near the origin for all times t (see eq.(23)) 
and thus the local maximum in po at x = 0 is preserved 
in p = Pol J . This non-generic behaviour occurs because 
of the lower smoothness of the spliced potential. The 
peak of p at X = 0 cannot be easily seen in the numerical 
solutions presented in the previous section because the 
height of the peak decreases exponentially in time. 

To identify the critical time Tc(p) one can apply the 
condition for a degenerate extremum, which is here 

J't'M/Mxo) = Po(xo)/po(xo). 

However, in the present case there is a more convenient 
approach, due to the fact that —J't/J increases mono- 
tonically in time for xq > 0. Because of the symmetry of 
the problem, we can restrict our attention only to posi¬ 
tive xo,x. To demonstrate the monotonicity for xq > 0, 
first note using the evolution equations (18) for J,K and 
similar evolution equations for J',K' [33] that 

^{K'J - J'K) = -V'''{x)J^. 
at 

Since particles in ID preserve their order, J > 0 and thus 
d/dt{K'J—J'K) < 0 for potentials of the sort considered 
here with V"'{x) > 0 for x > 0. Also, K'J — J'K = 0 at 
t = 0 so that 

^(-J'/J) = -{K'J - J'K)/> 0. 

Note that generally —Po/po > 0 at x > 0, for initial wave- 
packets with probability decaying monotonically away 
from the origin. Also, at time zero J'/J = 0, since 
Jo(xo) = 1. Thus, a simple procedure is to evolve —J'^/Jt 
in time and, as it rises upward from zero, identify the 
first time that its graph touches the graph of —pg/po. 

We can now apply these results to evaluate the split¬ 
ting time Tc{p) for the case of a Gaussian wave-packet. 
Since —Pq(xo)/po(xo) = xq/ct^, the condition for extrema 
of p becomes 

-J't{xo)/Jt{xo) = xo/cr^. 

It is most convenient to use this condition in a dimen¬ 
sionless form for uq = xo/i and r = t/{i^~°'/C)^^"^, or 

'^r(^o)/ dr('Uo) — p rto, 

with p = ^|<7. The full analysis of this equation for all p to 
find the mimimum r with solutions at positive uq requires 
a numerical study using the method discussed above [33], 
the results of which are presented in Fig.3. However, 
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10'® 10'^ 10'^ 10° 10° 10^ 10° 
/i 


(a) 



(b) 


FIG. 3 Splitting time Tc{^) plotted versus /i, (a) in 
log-linear and (b) in log-log, for both Kummer potential 
and Ist-order spliced potential. Notations as in Fig. 1. 


analytical results can be obtained asymptotically for /i ^ 
1 and /X 1. 

We first consider ^ ^ 1 for the Kummer potential. 
One can then just Taylor-expand to get 

Xt{xo) = xo + Xot + H- = Xo- ^V'{xo)t'^ H-■ 

We note here the convenient formulas [33] 


Differentiation using the formulas above and non- 
dimensionalization yields for r ^ 1 


J'riuo) 
•JT (tXo) 


1 

2 


W"(uo)t2 + • • •. 


Using further the Taylor expansion of the Kummer func¬ 
tion iUi(a, b]z) = l+^z-\ -, gives the result to linear 

order in and cubic order in uq. 

i(l - a)(l - = Q(1 - wo- 


Up to a rescaling of the variables, this is the standard nor¬ 
mal form of the pitchfork bifurcation, which is verified by 
the numerically obtained bifurcation diagram in Fig. 4. 
There is always the solution xxq = 0, but a symmetri¬ 
cal pair of nonzero solutions appears when the righthand 
side is positive, i.e. for t > Tc(/x) with 




((Tvlyi:) ' " ■« '■ 


This asymptotic behavior is verified in the numerical 
results plotted in Fig.3(b). Note that the linear scal¬ 
ing of Tc(/x) for small /x (but with a different prefactor) 
will hold for a generic, smooth symmetric potential with 
V"'{x) > 0 for a; > 0, by the same argument. We also ob¬ 
serve in Fig.3(b) the same linear scaling of Tc(/x) for the 
spliced potential, even though that potential has lower- 
order smoothness and the bifurcation in that case is an 
off-origin saddle-node type, as discussed in [33]. 



r 


u;(x) = -A^C^^-^x ■ iFi (^i^, 

VI'{x) = -A^cr-\F, i; 

U/"(x) = (1 - a)A„Ur-3x • iFi 

with the dimensionless constant 

A^ = 2(i+“)/°r((2 -b 0)12)1^/^. 


FIG. 4 Bifucation diagram for the Kummer potential with 
0 = 1/3 and for /x = 1, plotting uo for all extrema of the 
WKB position PDF versus dimensionless time r. Solid lines 
denote local maxima of the PDF, dashed lines local minima. 

We now turn to the limit /x 1. We have found no 
convenient analytic approach here for the Kummer po¬ 
tential, so that we consider instead the first-order spliced 
potential, where we can obtain closed expressions for Jt 
and J't/ Jt- For xq < ^ and t < ti{xo), one sees easily that 

C^IH 


Jt{xo) = cosh 


£{l-a)/2 


, J/(xo) = 0. 
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For xq < £ and t > ti{xo), we find that 

XoXt 


Jtixo) = 


-l+a^l_a 




- X, 


1+a ^ 


2C 
1 + a 


1/2 


{t - ti{xo)) 


Xo 


-l+a^l-a 


( ^ ) 

1/2 ^(l+a)/2y/y3^ 

\l + a) 



Finally, for xq > £ and all times t > 0 

1 


Jt{xo) = 


Xo 




Before proceeding, we note that these latter two expres¬ 
sions exhibit a jump-discontinuity in Jt(xo) at xq = £, 
which leads to the discontinuity in the WKB probabil¬ 
ity density p which was observed in the previous section. 
We also note here that Jt at long times appears to grow 
oc a;+ (t), but there is in fact a near-cancellation between 
two terms with such growth. The actual growth is 


Jt{xo) cx[x+{t)/£]*''^^°‘'^/‘^G{xol£), 

for some explicit function G(uo) [33]. This result has 
the interesting implication that the peaks in p(x^ t) at 
x±{t) have widths cx as t —)■ oo. However, 

the scaled density p{x,t) has peaks at ±1 with widths 
(X —)■ 0 as T —)■ oo. 

It is quite straightforward by differentiating the previ¬ 
ous results with respect to xq to obtain expressions for 
— J[{xo)/Jt{xo)- Complete formulas are given in [33], but 
here we write only the finite limiting result for r —>■ oo : 


*^00 (^o)/^oo ('^o) — T 


1 -t- Ha 


Uq 


(1 — a) -I- (1 -I- a)i 


Uo 


l-Wo + uly/l - ulF{uo) 


(27) 


with 


F{uo) 


/l 1-0 . 1 -i-3q . (1 - o) + (1-l-Q)«g 
1^2’ ~2(1-|-q)’ 2(1+ q)’ 2 


for Uo < 1, and 

-JUuo)/J^{uo) = ^^- ( 28 ) 

2 Uo 

for Uo > 1. This function is plotted in Fig.5. 

We note that the straight line p-^Uo for /i 1 inter¬ 
sects the curve —J'^{uo)/Joo{uo) at two non-zero values 
of Uo in the range 0 < uq < 1, one near uo = 0 and 
another near uo = 1. Because of the monotonic increase 
of — J'(uo)/Jt(uo) with T, there are no solutions at any 
time between these values of uq. We can determine the 
solutions near uq = 0 and uo = 1 by considering the 
behavior of — J^(uo)/Joo(uo) near these points. Near 
Uo = 0, 


^oo (llo)/ Joo{+0^ 


1 

Uo 



Uo 


FIG. 5 Plotted versus uo are the function — J(o(mo)/^oo(mo) 
for the spliced potential with a = 1/3, as solid line, and the 
linear function p^Uo for p = 5, as dashed line. 


SO the solution is at uq — 1/p and a height p. On the 
other hand, near uq = 1, 

c 


for a constant c, so the solution is at uo — 1 — c^/2/x^ 
and a height p^. Since the height of the first intersection 
is much smaller than the height of the second for /t 1, 
the graph of —J/(uo)/Jt(uo) touches the line p^uo for 
a smaller value of r near the first point uq — 1/p- This 
cannot happen at any time r < arccosh(l/uo), because 
— J/(uo)/Jr{uo) = 0 before then. This argument yields 
an exact low bound Tc{p) > arccosh(/x). In fact, the graph 
rises sufficiently rapidly after that time so that 

Tc{p) — arccosh(/t) ~ In/t, p+> 1] 

This asymptotic formula agrees with the numerical re¬ 
sults on splitting time Tc{p) for the spliced potential at 
p^ 1 observed in Fig.3(a). 

Furthermore, the same logarithmic asymptotics 
Tc{p) ~ In/i at large p appear in the plot for the Kummer 
potential in Fig.3(a). It is interesting that the same large- 
p asymptotics of Tc{p) hold for the pitchfork bifurcation 
in the Kummer potential and for the saddle-node bifurca¬ 
tion in the spliced potential. Although we have no proof 
of this asymptotics for the Kummer potential, there is a 
simple heuristic argument which leads to this conclusion 
in both cases. The regime pl^ 1 corresponds to an initial 
wavepacket with an extremely narrow spread in position 
space, a <+. £, very well localized in the region where the 
potential appears as an inverted oscillator. However, it 
is intuitively clear that the splitting of the wavefunction 
cannot occur until it has spread out of the parabolic re¬ 
gion of the potential. In the WKB limit the spreading 
is due to the classical trajectories, which escape that re¬ 
gion exponentially quickly. The time Tc{p) ~ In/i thus 


^oo 
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corresponds to the time it takes for classical particles at 
distance ~ a from the origin to reach a distance ~ 
These findings for rc(/i) are key results of our paper, 
providing the time of first onset of quantum spontaneous 
stochasticity phenomena. The implication is that in the 
WKB regime this is a short time of order 
The longest onset times are for /i ^ 1, but even here 
the growth of Tc(^) is only logarithmic, so that QSS phe¬ 
nomena are observable quickly unless the initial wave- 
packet spread is exponentially small, with cr = e~^I for 
M = In/i ^ 1. 


C. Scattering 


We now consider the problem of the scattering of a 
wavepacket with (x), (p) 0 when directed toward a 

UV-cutoff version of the cusp potential (7), which may 
be easier to study experimentally. To observe QSS in 
this setting it is clear that the classical dynamics for 
the initial condition ((x), (p)) must reach the phase-point 
(x,p) = (0,0), since only here is the flow vector field 
U(x,p) for the cusp potential non-Lipschitz. Unlike the 
problem of turbulent advection where the fluid velocity 
field is everywhere (nearly) non-Lipschitz, a very careful, 
fine-tuned choice of (x), (p) is here required to observe 
QSS. The only initial conditions for the unregulated cusp- 
potential which reach (0,0) are of the form (x+(t), V- (t)) 
and {x-{t),v+{t)) for some t > 0, where x±{t) are the 
extremal solutions given by (9) and v±{t) = dx±/dt. This 
condition imposes a relation between x and v: 

X =-sign(x) f j |x|^ (29) 

It is natural to expect that the pair ((x),(x)) must be 
chosen close to a point satisfying (29) in order to observe 
QSS also for a UV-cutoff version of the cusp-potential. 

To verify that this is so, we apply a simple criterion for 
presence of QSS. In order to have wave-packet splitting 
in a regulated cusp-potential it is enough to verify that 
the probabilities for the particle to be found to the left 
and right of the origin 

/ O f-\-oo 

dx \'ip{x,t)\‘^, p+{t)= dx |V’(a;,t)P, 

-oo J 0 


both remain strictly positive for all times t > ts,, the 
splitting time, in the limit £, cr —>■ 0. One expects correc¬ 
tions to of order 0((£^““/C')^/^. These probabilities 
are easy to evaluate in the WKB limit as 


P-ii) 




P+{t) 




assuming for simplicity a minimum-uncertainty Gaussian 
initial wavepacket, where 


y*{i) 


,i.„ ( "<‘(0)-M 

y cr 


(30) 


and x)"^(0) = xq for the initial condition (xq, {v)) which 
arrives to the origin 0 under the classical dynamics ex¬ 
actly at time t. It follows that, to see wavepacket splitting 
within WKB, one must have \y*{t)\ < oo for all t > t*. 

We can calculate the above quantities analytically for 
the spliced potential. Note that x)"^(0) must lie between 
the origin 0 and x^o = xJ^^(O), the point at which a 
particle with initial velocity (v) arrives to the origin only 
after an infinite time. Energy conservation gives Xoo as 


1 + a 2 2 


l+a 


1 — a 
1 + a 


since the particle must arrive to 0 with vanishing veloc¬ 
ity. If we introduce the position x* which satisfies the 
constraint (29) with the velocity (v), then this becomes 

ia;oor+“ = ix*r+“ + i(i-a)£i+“. (31) 

As expected, Xoo —>■ x* as £ —)■ 0. The above derivation 
has assumed that |(x)| > so that |xoo| > £■ If 

instead, |(x)| < then the point Xoo lies in the 

inverted-oscillator region of the spliced potential. In that 
case it is easy to show again by energy conservation that 


Xoo = -{v)hi 
with 7 = (C'/£i-“)i/2. 

The classical solution which starts at Xoo must spend 
an infinite amount of time in the inverted-oscillator re¬ 
gion |x| < £ before reaching 0. This can be verified by 
exactly solving the dynamics in the inverted-oscillator 
region as 


X = xo cosh( 7 <) H —Vo sinh( 7 t), 

7 

so that the particle reaches the origin when |xo| > 7 |xo| 
and xo,xo have opposite signs, with this process taking 
a time t = (l/7)arctanh|7Xo/wo|- To apply this result to 
solutions starting at Xq with |xo| > we note that such 
solutions reach distance £ from the origin with a velocity 
ve(xo) given by energy conservation as 


2C 

1 -I- a 


Ja;* 


l+a 


1^0+“ , 


or, using (31), as 


= {.l£f + - |a:or+“] • 

Hence the time that classical solutions for such initial 
positions xq spend in the inverted-oscillator region is 

X / 'y£ 

Teixo) = — arctanh —-—r 

7 V^ua^o) 

For the special case Xq = x^o, one has V( = ^£ and thus 
the time is infinite, as expected. 
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T = 0 



(a) 


r = 10.00 



(c) 


T = 5.00 



(b) 


T = 15.00 



(d) 


FIG. 6 WKB position-space probability densities for an incoming wave-packet with {v )/= 20, {x) = x^o, and 
a = ^ at dimensionless times r = 0, 5, 10, 15. Here a = 1/3 and notations are as in Fig. 1. For movie, see [33]. 


(/x^{t) = 4.4 X 10“’ = 4.35 X 10“® 




i/x,{t) = 4.35 X 10-^ e/x+{t) = 4.35 x 10-® 




(c) 


(d) 


FIG. 7 WKB position-space probability densities for incoming wave-packet with (v), {x} and time t as in Fig. 6(d), but now 
with i reduced to give ilxc{t) = 4.35 x 10“^, 10“®, 10“^, 10“®, with notations as in Fig. 1. For movie, see [33]. 


These considerations fully characterize the value xq = (0). The important constraint comes from the solu- 
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tion ( 11 ) for the dynamics in the cusp-region |x| > I 
of the potential. Using this equation for final position 
\x\= £ and time t — ti{xq) to reach that point yields 


^(i-a)/2^^ / __(!-«) . l + 3a 


Xq 


2’ 2(1-h a) ’ 2(1-ha) ^ 


2’ 2(1-ha)’2(1-ha)’ 


1 + Q: 

a^o 


l-l-aN 


= - 2 (l-«) 


2 C 
1 -h a 


1/2 


{t-Tt{xt))). (32) 


Here we have assumed a negative-energy solution appro¬ 
priate to a point near Xco and taken a minus sign on the 
righthand side for the case sign(a:o-(u)) = —I- To evaluate 
xo = a;)"^( 0 ) asymptotically for small i it is convenient to 
parameterize this solution by the dimensionless quantity 
^ > 0 introduced through the formula 


Fo 


= \Xc 


-it 




In 


terms of this parameter vf = C and 


Ti{i) = — arctanh 


1 -h a 

l-ha-h2C' 


(34) 


Also the quantity xq appearing in (32) must be chosen 
according to the condition 


which can only hold if t—re = t*-hO((£^ “/C)^/^), where 
is the time given by the relation 


,(l-a) 


2C 
1 -h a 


1/2 


2/(l-a) 


or the time t* required for the extremal solutions (9) of 
the cusp-potential to transit from a;» to 0. We conclude 
that Xq = x^^{0) for the spliced-potential is given by 

|xor+“ = |xoor+“ - 2(1 + a)e-2'^(‘-**)£i+“ (36) 

in the limit f —> 0, for Xoo in (31) and 7 = (Cj■> 
for all t > The simple picture that emerges here is that 
solutions starting at xq = a;)"^(0) given by (36) reach a 
distance £ from the origin in a time ~ t*, nearly indepen¬ 
dent of £, and then take an additional time t = t — ts, 
to reach the origin, for a total time t. This is made pos¬ 
sible by the fact that extremal solutions (9) of the cusp- 
potential transit between x* and 0 in a finite time t*, as 
a consequence of the non-smoothness of that potential. 

With these results in hand, we can now discuss the 
wavepacket scattering in the WKB limit. With the mean 
velocity (v) fixed, there are various possible choices of 
(x) to observe splitting of the wavepacket and QSS. For 
example, one may choose (x) = Xoo- Since (36) implies 


i?n = - 


C 


|l+a _ 1„2 _ ^ ( 1+0 


T—kol T 

1 -b a 2 1 


+ a 


which gives 


ixor+“ = 




fl+0 


(35) 


With these choices xq = ^^^^(O) is the solution of the 
fixed-point condition (32). 

We can now evaluate the solution asymptotically as 
£ —> 0. Because 7 —> 00 in that limit, the only way to 
satisfy (34) for a fixed Hnite value Tf(^) = r > 0 is to 
choose the parameter i exponentially small: 

^ = 2(1 + a)e-2'^U 




0 

1 + Q 

^00 

1 - 

X:^ 



almost any reasonable choice of ct, e.g. a (x £^ for any 
power p, will give 


y*{t) = 0, t > t». 

This corresponds to a symmetrical splitting, with both 
P- (t) = p+ (t) = I for t > . Another natural choice is 

(x) = X*. Using (36),(31), it follows that 


x^ ^(0) = X* 


^ 1 — a 

£ 

1 + Q 

1 + / 



2(1 -b a) 

X* 



This implies that xq in (33) must be exponentially close 
to Xoo and, in particular, xq —t x*, verifying our original 
expectation. Furthermore from (35) 


£ 


1 + Q: 




a 


and from (35),(33) 



1+a / 

i 

Xo 

= o( 

Xo 

V 

X:^ 



^ 0 . 


These results in the fixed-point condition (32) give 


so that taking a = +^+“/|x*|“ gives 
/ \ 1 1 — a . , , 

y»(t =^——sign(x* , t>U. 

2 p 1 -b a 

Hence p+{t) = |erfc( ^;(^ 2 ^ l^) t in terms of 
the complementary error function, corresponding to an 
asymmetrical splitting. Still another possible choice is 
(x) = X* — KCT for any cr vanishing more slowly than 
so that —>■ 0 as £ —>■ 0. In this case 

y*{t) = K, t > t*. 

and p+{t) = +rfc(«;) for t > This would hold with 
the choice a = £!p oi the previous section, for p fixed. 
Although there seem to be a plethora of possibilities, they 
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are actually quite restricted since one must always make 
a careful fine-tuning so that (x) —)■ cc* as f > 0 in order 
to observe QSS in the scattering process for fixed {v). 

We now present numerical WKB results for scatter¬ 
ing of a wave-packet off both the Kummer potential and 
the spliced potential. Here we shall consider only the 
choice (x) = x^o for fixed (v) and a = £, based upon 
the analytical results for the sliced potential. Using the 
same choices for the Kummer potential is a good test on 
the robustness of the derived selection criteria for small 
changes in the regularization of the potential. (Note that 
one could also make an alternative choice of Xoo for the 
Kummer potential by using its own energy conservation 
law to choose the position of the particle with velocity (v) 
which arrives to the origin with zero velocity.) We con¬ 
sider first the choice a = 1/3 and = 20 so 

that Xoo/£ = —66.05 and r* = = 9.902. 

In Fig. 6 we show the time-evolution of the WKB posi¬ 
tion PDF by plotting it at four successive choices of the 
dimensionless time r. The classical solution with which 
we now compare is Xc(t) = a;±(|t —t*|), ± = sign(t —t*), 
which is incoming as X- and outgoing as x+. We see that 
the PDF’s for both potentials evolve very similarly and 
split at a time near . It is notable that in the final frame 
for r = 15, the probabilities p-,p+ are nearly equal to 
the symmetrical value 1/2 but the shapes of the peaks 
on the left and right are different. Next we study the 
effect of decreasing values of cr = ^. In Fig. 7 we fix {v), 
(x) = Xoo, final time t, and position Xc{t) from the last 
panel in Fig. 6 and successively reduce £/xc{t). There is 
clear evidence of QSS, with the PDF converging to a pair 
of delta-functions. Even though p+,p- —)■ 1/2, the shape 
asymmetry of the PDF peaks remains. The most impor¬ 
tant observation is that the results for the Kummer and 
spliced potentials are quite similar, giving hope that our 
selection criterion, although derived rigorously only for 
the spliced potential, will apply more generally. Other 
choices of the fine-tuning of the initial wave-packet are 
implemented by the codes in [33]. 


V. NUMERICAL SCHRODINGER EVOLUTION 

Thus far, we have considered a WKB regime in which, 
formally, the limit h/m —>■ 0 is taken first with position 
spread Axq = cr of the wavepacket fixed and then subse¬ 
quently letting cr —)■ 0. However, this is not the only limit 
in which deterministic classical dynamics is generally ex¬ 
pected to emerge. One can also let h/m and Axg go to 
zero together, e.g. with position spread Axg oc ih/m)^ 
and velocity spread Auq oc (/i/m)^“^ for any 0 < /3 < 1, 
consistent with the Heisenberg uncertainty relations. For 
general smooth potentials and for an initial minimum- 
uncertainty wavepacket with Axq, Auq chosen as above, 
it can be shown that the position and velocity spreads 
Axt, Avt at later times t vanish for the limit h/m —)■ 0 
and the mean values {x)t, {v)t evolve according to the 
classical dynamics [46]. Note that WKB corresponds to 


the special case j3 = 0. However, if the potential is rough 
down to a length-scale £ and if £ —)■ 0 as h/m —>■ 0, then 
particle motion could remain stochastic in any of these 
alternate classical limits as well. 

To give a careful formulation of such alternative semi- 
classical limits, we define a “dimensionless-?!” parameter 

^ ^ TOC'V2^(3-ea)/2 ’ 

where 6 is some suitable macro-length scale, which could 
be the outer length L of the potential, the length D of 
the space domain, or some specified fraction of either of 
these lengths. The quantity defined in (37) thus corre¬ 
sponds to the “inverse Peclet number” k/UL in terms 
of the turbulent advection - quantum mechanics analogy 
discussed in Section HI. The alternative semi-classical 
limits are then defined by taking e <C 1 with 

a = £/pL, £ = 6 ■ e^, 0 < /3 < 1 . 

This choice of £ gives the potential additional mass- 
dependence beyond that in the prefactor C ex. 1/m. If 
the Schrodinger equation (19) is non-dimensionalized by 
introducing x = x/5, i = t/(5^““/C')^/^, then it becomes 

iedi't/=+ V{x)t/ (38) 

with now m = C = 5 = 1 and £ = for e ^ I. 

We numerically study the semi-classical limit in the 
setting described above with (3=1/2 and 6 = D/12?>'k. 
We modified a publicly available code [47] to solve the 
Schrodinger equation in one dimension for the Kummer 
potential (12), in the dimensionless formulation (38). 
The code pySchrodinger employs the so-called Strang 
operator-splitting spectral method, which is especially 
suitable for the semi-classical limit problem [48]. This ap¬ 
proximation introduces second-order time-discretization 
errors by time-splitting, but has exponential accuracy 
in spectral computation of spatial derivatives for suffi¬ 
ciently smooth wave-functions. It has the great advan¬ 
tage of giving evolution which is unitary, time-reversible, 
and gauge-invariant under spatially constant shifts in the 
potential, to machine precision. In the numerical ex¬ 
periments presented here we study an initial minimum- 
uncertainty wave-packet with {x) = (u) = 0 and p = 1, 
making a sequence of runs with e = 2“" for all integers 
—7 < n < 10. We resolve both the spatial and temporal 
grid to 0(e) and performed refinement studies to verify 
that our numerical simulations are converged. 

In Fig. 8 we show for e = 2“^ a sequence of snapshots 
at successive times of the position-space PDF, along with 
the positions of the extremal classical solutions. Just as 
in the WKB limit, the PDF splits into a bimodal distri¬ 
bution at an early time and with the two peaks closely 
tracking the extremal classical solutions. 

We can also study splitting in momentum-space by 
taking the Fourier transform of the wave-function, de- 
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FIG. 8 Position-space PDF for numerical Schrodinger solution with e = 2 at successive times i — 0, 2, 4, 6. Asterisks on 
the abscissa denote positions of the extremal classical solutions. For movie, see [33]. 
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FIG. 9 Momentum-space PDF for numerical Schrodinger solution with e = 2 at successive times t = 0, 2, 4, 6. Asterisks 
on the abscissa denote momenta of the extremal classical solutions. For movie, see [33]. 


fined in dimensionless form by 


-tpipj) 




and implemented numerically by a fast Fourier transform 
(FFT). The momentum-space PDF p{p,i) = is 

plotted in Fig. 9 for the same cases as in the preced¬ 
ing Fig. 8. For comparison, we also plot the momenta 


\p2/K^ 
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e = 2~' e = 2“° 



£ = 2“* e = 2"“ 



(c) (d) 


FIG. 10 Position-space PDF for Schrodinger solution at t = 6, for e = 2 ^,2 ®, 2 2 with notations as in Fig. 8. For 

movie, see [33]. 
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FIG. 11 Momentum-space PDF for Schrodinger solution at t = 6, for e = 2 ^ ,2 ®, 2 2 with notations as in Fig. 9. 

For movie, see [33]. 


p±{t) = mv±{t) of the extremal classical solutions in 
non-dimensionalized form. We see that there is split¬ 
ting of the wave-packet also in momentum-space and the 
two peaks of the momentum-space PDF closely track the 


momenta of the extremal classical solutions. Although 
we did not show it earlier, one can obtain similar results 
for the WKB limit. Numerically, the WKB wavefunction 
can be interpolated to a regular space grid with spacing 
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0{e) and then transformed by FFT. One observes the 
same splitting in momentum-space there. 

We also study the behavior of the position-space and 
momentum-space probability distributions at a fixed 
time t and £ = a (x for a decreasing sequence of 
£ values. The results are presented in Figs. 10-11. In the 
semiclassical limit e ^ 1, the probability densities in¬ 
creasing concentrate around the positions and momenta 
of the two extremal classical solutions, providing good 
evidence of QSS in both position-space and momentum- 
space. We can quantify this further by calculating the 
position- and momentum-space uncertainties as a func¬ 
tion of e at the same fixed time t. These results are shown 
in Fig. 12. The horizontal lines are the corresponding 
theoretical predictions from (26) for the classical dis¬ 
persion. We see a definite tendency of the numerical 
uncertainties to converge to their theoretical limit and, 
most importantly, become asymptotically independent of 
e, corroborating the visual evidence in Figs. 10-11. 



c 

FIG. 12 Uncertainties of position and momentum at the 
fixed time t in Figs. 10-11, as functions of e. The heavy solid 
curve is the numerical position-space dispersion as defined 
by Eq. (25) and the light solid line is the position dispersion 
of the two classical extremal solutions (26). The thick 
dashed curve is the numerical momentum-space dispersion 
and the thin dashed line is the momentum dispersion of the 
two classical extremal solutions. 

The previous results show that quantum systems may 
remain stochastic (“God still throws dice”) in the clas¬ 
sical limit, where determinism is expected. However, 
it might appear, even more shockingly, that the par¬ 
ticle models discussed here become increasingly accu¬ 
rately described as two-state quantum systems in the 
limit h/m —> 0, with just two possible values for the 
position (or momentum). If so, then it could be possible 
by experimental design of such systems to construct a 
massive “Schrodinger cat” with arbitrarily large ml The 
key question here is whether it is possible to achieve a 
coherent superposition of the two different states of the 
particle, leading to quantum interference effects. If there 
is no such coherence, then any linear superposition of the 
two possible states really corresponds just to a statistical 


mixture within classical probability theory. 

There are compelling reasons for lack of coherent su¬ 
perpositions, arising from the quantum uncertainty rela¬ 
tions. Note that coherence would require separately for 
each branch of a “split” wavefunction that the position- 
space spread be comparable to the de Broglie wavelength. 
Ax ~ A. Otherwise, the phase of each branch wavefunc¬ 
tion would oscillate rapidly within its envelope of width 
Ax and the two branches, even if spatially recombined, 
would have overlap nearly zero. In that case, no inter¬ 
ference between them could ever be observed. However, 
choosing Ax ~ A oc h/m in each branch leads by the 
uncertainty relations to Av = 0(1) for each branch in¬ 
dividually and thus very rapid spatial spreading. In the 
limits where classical dynamics is expected. Ax ^ A, and 
thus no coherence of the branches. Hence, the QSS phe¬ 
nomenon allows indeterminism in the classical limit but 
not quantum superpositions or entanglement. 

These arguments can be verified using our numeri¬ 
cal Schrodinger solutions. Although the two individual 
peaks of the position-space PDF in Fig. 10 become in¬ 
creasingly narrow oc as £ —> 0, they in fact con¬ 
tain increasingly rapid oscillations on the scale of the de 
Broglie wavelength ~ 0{e). We show this in Fig. ?? by 
zooming in on one of the two peaks (the rightmost) and 
plotting the function |Re('0)p = pcos‘^{(p), whose enve¬ 
lope is the position-space PDF but which oscillates with 
the complex phase Lp of the wave-function. Clearly, the 
wave-functions become extremely oscillatory within the 
envelope of a single peak as e is decreased. This rapid 
phase oscillation renders the two branches of the wave- 
function decoherent, as it becomes impossible in practice 
to recombine them spatially in phase. This is also a fea¬ 
ture of the WKB wave-function studied previously, as 
can be seen by calculating its phase St{x), the classical 
action (16), and observing that it grows with increasing 
mass m of the particle. 

VI. POSSIBLE EXPERIMENTS 

Such indeterministic, stochastic behavior persisting in 
the standard classical limits should be observable in lab¬ 
oratory experiments. We have examined in this paper 
only the simplest ID model which exhibits the QSS phe¬ 
nomenon and it may not be the easiest to situation to 
study experimentally. However, the results will readily 
extend to higher-dimensional cases that may be easier to 
create in the laboratory. The main requirement is that 
the repulsive potential of the quantum particle should be 
(nearly) non-smooth at the maximum energy location. 
For example, the radially symmetric version of the cusp 
function (7) as potential energy per unit mass, i.e. 

V(x) =r = |x|, 

1 + a 

will lead to the same QSS phenomena. By standard argu¬ 
ments, the 3D Schrodinger equation for an angular mo- 






19 





15 20 25 30 35 

X 


(d) 


FIG. 13 Plot of |Re(^)P near the right peak of the position-space PDF, for e = 2^, 2^, 2“^, 2“^, as light black line. The 
solid black curve is the plot of the position-space PDF p, the envelope of the first graph. Each plot window height is a fixed 
multiple of the maximum value of the peak and each window width is a hxed multiple of the peak half-width. 


mentum £ initial wave-packet will reduce to a ID equation 
with an effective potential 14ff(r) = V{r) -\- 
particular, for an initial £ = 0 wavepacket concentrated 
at the origin, the analysis in this paper carries over and 
predicts a spherical wave of probability amplitude ra¬ 
diating from the origin in the standard classical limits 
where one naively expects the particle to remain sitting 
at the unstable maximum at the origin. In this case, 
the randomness involves a continuous infinity of random 
outcomes, corresponding to the possible directions of the 
outgoing particle, not just two as in the ID model. 

To verify the phenomenon as QSS, there are key signa¬ 
tures for the experimentalist to check. First, there should 
be super-ballistic power-law spreading of the wavepacket 
as in (26). Second, this spreading should be indepen¬ 
dent of the initial width of the wavepacket a and in¬ 
ner cutoff £ of the potential as those are decreased (but 
still much larger than the de Broglie wavelength). Third, 
the spreading should be independent of h/m, which may 
be testable in some experiments by considering parti¬ 
cles with differing masses. (Actually in a fixed potential, 
there will be a particular mass-dependence, through C or 
classical inertial effects, rather than strict independence.) 
It may be easier for experimentalists to observe the cor¬ 
responding signatures for momentum-space spreading, as 
discussed in section V. Any experiment attempting to 
observe QSS will need low temperatures and low levels 
of other sources of noise, to ensure that quantum fluctu¬ 
ations are driving the phenomenon. We consider here 


briefly three possible experimental realizations: ultra¬ 
cold AMO, ultra-cold neutrons, and optical analogues. 

Ultra-cold atomic-molecular-optical (AMO) systems 
exhibit quantum evolution of neutral particles in an 
electric dipole potential created by the AC Stark effect 
[49, 50]. This potential has the form 


U(x) 


Sttc^ r 

2u;§ A 


/(x) 


where cuq is the optical transition frequency, F is the 
damping rate corresponding to the spontaneous decay of 
the excited level, I (x) is the spatial intensity of the laser 
light with frequency w, and A = w — wo is the detuning 
parameter. Achievable potentials have only a small en¬ 
ergy scale (mK), but there is great flexibility in designing 
their spatial form (e.g. see [51, 52]). Attractive electric 
dipole potentials are widely used to create particle traps 
for atomic bose-condensates and for quantum computing. 
A stationary repulsive potential for study of QSS could 
be accomplished with a blue-detuned (A > 0) focused 
beam or red-detuned (A < 0) hollow beam. The QSS 
phenomenon requires a large scale ratio, L/£ ^ 1, say, 
L/£'^ 100, which could be achieved perhaps with large- 
aperture lasers or a large number of individual beams. 
To get the cusp-potential would require a fine control of 
the spatial intensity distribution /(x), so that a power- 
law such as /(x) oc |x|^+“ could be obtained to a suit¬ 
able approximation. The effects described here should 
not require an exact power-law but could have discrete 
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steps of intensity or some rapid variation that averages 
out on the dynamical times scales of the particle. To 
achieve such a spatial distribution, one might use meth¬ 
ods employed for quantum computing applications, such 
as artificial holograms with acousto-optic deflectors [51], 
spatial light modulators [52], or other devices. In con¬ 
trast to quantum computing applications, one does not 
require rapid manipulation and control of potentials, but 
instead accurate specification of the intensity spatial dis¬ 
tribution over a broad enough region. This seems a very 
promising approach to explore QSS in the laboratory. 

Another possible system with which to realize QSS is 
ultra-cold neutrons [53-55]. These have energies 1-1000 
neV, temperatures ~mK, de Broglie wavelengths A = 
10-1000 nm, velocities n ~ 5 m/s. Neutrons interact 
through through all four fundamental forces —strong, 
weak, electromagnetic and gravitational—which provide 
various methods to generate potentials acting upon them. 

One possibility is the magnetic dipole potential due 
to an external magnetic field, l^(x) = —/x*B(x), which 
acts on the neutron dipole moment with /x = 60 neV/T. 
(Magnetic dipole forces are also an option for AMO sys¬ 
tems). Magnetic field strengths B as large as 30 T are 
currently achievable with superconducting magnets and 
generally somewhat smaller strengths ~ 1 T for per¬ 
manent magnets. Unless magnetic fields change very 
rapidly, the alignment of neutron polarization generally 
follows adiabatically the local magnetic direction, with 
no spin flips. Thus the dipole interaction is repulsive 
for spins aligned or magnetic moments anti-aligned with 
the field (low-field seeking, LFS, state), and attractive 
for spins anti-aligned/moment aligned (high-field seek¬ 
ing, HFS, state). Magnetic bottles for spin-polarized neu¬ 
trons, as originally proposed in [56], have been achieved 
for both LFS and HFS spin states and for electromagnetic 
and permanent magnets. See e.g. [57] for many refer¬ 
ences. One could, in principle, use similar methods to 
generate a (UV-regulated) cusp potential. An advantage 
of this approach is that magnetic interactions are ideal 
conservative and relatively “clean” compared to methods 
involving material mirrors/walls. However, it is difficult 
for us to conceive of a practical arrangement of perma¬ 
nent magnets and electric currents that would produce a 
cusp-like potential, with sufficient accuracy. 

It is also possible to exploit gravity, since ultra-cold 
neutrons have kinetic energies comparable to their grav¬ 
itational potential energies (102 neV per meter height). 
This can be done in conjunction with material mirrors 
(essentially, the strong interaction with the atomic nuclei, 
modelled by nuclear optical potentials). Such a gravity 
-I- strong force approach has already been used to create 
attractive potential wells and quantum bound states for 
ultra-cold neutrons [58]. For classical dynamics the ID 
cusp potential can be generated by a particle constrained 
to move along a suitable space curve (“frictionless bead 
on a wire”) falling under gravity, so long as the required 
forces are < mg. Quantum mechanically one might be 
able to exploit micron-scale or submicron-scale channels 


as neutron guides. Polycapillary glass fibers consisting of 
bundles of ~ 1000 individual lead-silica glass capillaries, 
each with diameter ~ 6/xm, have been demonstrated to 
act as guides for thermal neutrons [59] and are commer¬ 
cially available [60]. It is conjectured that neutron guides 
with inner diameters as small as 5-7 nm may be possible, 
e.g. using carbon nanotubes [61, 62]. When the diame¬ 
ter of the guide is comparable to the wavelength of the 
confined ultra-cold neutron, then simple energetic con¬ 
siderations as well as more detailed analysis [63, 64] im¬ 
ply that the wave propagation should be quasi-lD, with 
little excitation of high-order modes transverse to the 
tube axis. Using gently bent guide tubes to reduce losses 
[59, 63] and considering times short enough so that neu¬ 
trons have velocities sufficiently small to remain confined, 
the ID models of this paper might be obtained. 

Another entirely different class of experiments that 
could realize our ID models are based on optical ana¬ 
logues of Schrodinger evolution [65, 66]. For example, 
spatial propagation of a monochromatic light beam in 
guiding dielectric structures which are weakly curved 
along the propagation direction are described by the 
paraxial optical wave equation: 

Here "0 is the complex electric field amplitude modified 
by reference-frame and gauge transformations, z is the 
paraxial propagation distance, % = \/2tt is the reduced 
wavelength, V(x) = [ng—n?{x)]/{2ns) — ns — n(x), n{x) 
is the refractive index profile of the guiding structure, 
is the reference (substrate) refractive index, xo{z) is the 
axis bending profile, and F{z) = —n‘^Xo{z). For full de¬ 
tails, see [65]. It appears feasible to produce a repulsive, 
cusp-like potential by manufacturing a dielectric medium 
with a suitable specification of n{x) < Ug. Such experi¬ 
ments are not truly quantum-mechanical, however, but 
are to be more properly considered as “analogue simu¬ 
lations” of the Schrodinger equation. Their value is to 
be assessed by comparison with direct numerical simu¬ 
lations of the Schrodinger equation, like those presented 
in this paper. Optical analogue experiments are useful if 
they can probe a complementary (or wider) range of pa¬ 
rameters than those achievable by computer simulations. 

The specific suggestions made here for experimental re¬ 
alization of QSS may not be the most feasible ones and, 
in any case, are mere sketches of possible approaches. 
Other possibilities will surely occur to our ingenious ex¬ 
perimentalist colleagues. However, we believe that our 
theoretical analysis makes a compelling case that such 
breakdown of classical determinism will eventually be ob¬ 
served for quantum systems in the laboratory. For this 
not to be true, quantum mechanics itself would need to 
fail in the distinctive situations considered here. 



21 


VII. CONCLUSIONS 

This paper has developed a novel prediction of quan¬ 
tum mechanics: the breakdown of determinism in the 
standard classical limit (large mass, small initial spreads 
in position and velocity) when the force-field acting on 
the quantum particle is nearly rough. Like the classical 
spontaneous stochasticity effect, it is produced by an ex¬ 
plosive, super-ballistic acceleration of the particle and the 
resultant “forgetting” of its initial conditions. Quantum 
spontaneous stochasticity is, in some sense, a diametric 
opposite of Anderson localization [67]. That effect sup¬ 
presses spreading of the wave-function in a random (dis¬ 
ordered) potential, whereas QSS corresponds to an ac¬ 
celerated spreading. Furthermore, Anderson localization 
is essentially a wave-phenomenon, arising from destruc¬ 
tive interference between multiple-scattering paths. QSS 
is instead an essentially classical particle phenomenon: 
although the origin of the randomness lies in the quan¬ 
tum fluctuations, quantum stochasticity is magnified to 
macroscopic scales by classical non-smooth dynamics. 

In this paper we have confined our attention to the 
initial-value problem for a quantum particle, described 
by the non-relativistic Schrodinger equation, but QSS 
should appear in other guises and theoretical frameworks. 


The Wigner distribution function method, previously ex¬ 
ploited in [28], appears to be a powerful method to an¬ 
alyze simultaneously the position-space and momentum- 
space spreading of the quantum particle. There should 
also be very interesting manifestations of classical non¬ 
smooth dynamics in the Feynman path-integral represen¬ 
tation of quantum scattering amplitudes [27], not only 
in the classical limit but also in the fully quantum- 
mechanical setting. Studying the quantum phenomenon 
will help to shed further light on classical spontaneous 
stochasticity, which has fundamental importance for as¬ 
trophysics, geophysics, and even engineering applica¬ 
tions. It is also critical to continue testing quantum me¬ 
chanics in novel regimes that were previously unexplored. 

ACKNOWLEDGMENTS 

We would like to thank Eberhard Bodenschatz, Kurt 
Eyink, Igor Kolokolov, Jorge Kurchan, Guy Marcus and 
Thomas O’Conner for very useful conversations and sug¬ 
gestions on possible experimental realizations of quantum 
spontaneous stochasticity. Abraham Harte informed us 
of the philosophical paper by Norton. TD would like to 
thank Cristian C. Lalescu and Navid Constantinou for 
their help with the Schrodinger simulations. 


[1] P. Hartman, Ordinary Differential Equations: 2nd Ed., 
Classics in Applied Mathematics (Society for Industrial 
and Applied Mathematics, SIAM, 3600 Market Street, 
Floor 6, Philadelphia, PA 19104, 1982). 

[2] J. Norton, Philosopher’s Imprint 3, 1 (2003). 

[3] D. Bernard, K. Gaw§dzki, and A. Kupiainen, J. Stat. 
Phys. 90, 519 (1998). 

[4] M. Chaves, K. Gaw§dzki, P. Horvai, A. Kupiainen, and 
M. Vergassola, J. Stat. Phys. 113, 643 (2003). 

[5] L. F. Richardson, Proc. R. Soc. Lon. Ser-A 110, 709 
(1926). 

[6] A. Lazarian and E. T. Vishniac, Astrophys. J. 517, 700 
(1999). 

[7] E. Vishniac and A. Lazarian, in Plasma Turbulence 
and Energetic Particles in Astrophysics; Proceedings 
of the International Conference, Cracow, Poland, 5- 
10 September, 1999, edited by M. Ostrowski and 
R. Schlickeiser (Obserwatorium Astronomiczne, Uniwer- 
sytet Jagiellohski, Krakow, 1999) pp. 182-189. 

[8] G. L. Eyink, A. L. Lazarian, and E. T. Vishniac, Astro¬ 
phys. J. 743, 51 (2011). 

[9] G. Eyink, E. Vishniac, C. Lalescu, H. Aluie, K. Kanov, 
K. Biirger, R. Burns, C. Meneveau, and A. Szalay, Na¬ 
ture 497, 466 (2013). 

[10] G. L. Eyink, C. R. Phys. 7, 449 (2006). 

[11] G. L. Eyink, Phys. Rev. E 83, 056405 (2011). 

[12] R. Bitane, H. Homann, and J. Bee, J. Turbul. 14, 23 
(2013). 

[13] M. Bourgoin, N. T. Ouellette, H. Xu, J. Berg, and E. Bo¬ 
denschatz, Science 311, 835 (2006). 

[14] A. N. Kolmogorov, Dokl. Akad. Nauk SSSR 30, 299 


(1941). 

[15] L. Onsager, Nuovo Cimento Suppl. 6, 279 (1949). 

[16] G. L. Eyink and K. R. Sreenivasan, Rev. Mod. Phys. 78, 
87 (2006). 

[17] B. 1. Shraiman and E. D. Siggia, Phys. Rev. E 49, 2912 
(1994). 

[18] G. L. Eyink and T. D. Drivas, J. Stat. Phys. 158, 386 
(2015). 

[19] R. H. Kraichnan, Phys. Fluids 11, 945 (1968). 

[20] G. Falkovich, K. Gaw§dzki, and M. Vergassola, Rev. 
Mod. Phys. 73, 913 (2001). 

[21] K. Gaw§dzki and M. Vergassola, Physica D 138, 63 

( 2000 ). 

[22] W. E and E. Vanden-Eijnden, Proc. Natl. Acad. Sci. 97, 
8200 (2000). 

[23] W. E and E. Vanden-Eijnden, Physica D 152-153, 636 

( 2001 ). 

[24] Y. L. Jan and O. Raimond, Ann. Prob. 30, 826 (2002). 

[25] Y. L. Jan and O. Raimond, Ann. Prob. 32, 1247 (2004). 

[26] T. D. Drivas and G. L. Eyink, J. Fluid Mech. , to be 
submitted. 

[27] R. P. Feynman, Rev. Mod. Phys. 20, 367 (1948). 

[28] A. Athanassoulis and T. Paul, Mathematical Models and 
Methods in Applied Sciences 22, 1250038 (2012). 

[29] T. Paul, in Microlocal Methods in Mathematical Physics 
and Global Analysis (Springer, 2013) pp. 49-52. 

[30] R. Bliimel, Acta Phys. Pol. A 93, 7 (1998). 

[31] L. Couchman, E. Ott, and T. M. Antonsen Jr, Phys. 
Rev. A 46, 6193 (1992). 

[32] R. Bliimel, T. Antonsen Jr, B. Georgeot, E. Ott, and 
R. Prange, Phys. Rev. E 53, 3284 (1996). 



22 


[33] See Supplemental Materials at, littp://link.aps.org/ sup¬ 
plemental/??. ????/PhysRevX. ?.?????? 

[34] M. Abramowitz and I. Stegun, Handbook of Mathemati¬ 
cal Functions: with Formulas, Graphs, and Mathemati¬ 
cal Tables, Dover Books on Mathematics (Dover Publi¬ 
cations, 2012). 

[35] Bateman Manuscript Project, A. Erdelyi, and H. Bate¬ 
man, Higher Transcendental Functions, vol. 1 (Dover 
Publications, 2007). 

[36] Bateman Manuscript Project, H. Bateman, and 
A. Erdelyi, Tables of Integral Transforms: Based, in 
Part, on Notes Left by Harry Bateman, v. 1 (McGraw- 
Hill, 1954). 

[37] R. H. Kraichnan, Adv. Math. 16, 305 (1975). 

[38] A. Lazarian, E. T. Vishniac, and J. Cho, Astrophys. J. 
603, 180 (2004). 

[39] D. Ryu, H. Kang, J. Cho, and S. Das, Science 320, 909 
(2008). 

[40] A. Lazarian, Space. Sci. Rev. 181, 1 (2014). 

[41] K. Gottfried and T. Yan, Quantum Mechanics: Fun¬ 
damentals, Graduate Texts in Gontemporary Physics 
(Springer, 2003). 

[42] L. O’Raifeartaigh and A. Wipf, Found. Phys. 18, 307 
(1988). 

[43] Y. Goldfarb, J. Schiff, and D. J. Tannor, J. Chem. Phys. 
128, 164114 (2008). 

[44] G. Barton, Ann. Phys. 166, 322 (1986). 

[45] G. Yuce, A. Kilic, and A. Coruh, Phys. Scripta 74, 114 
(2006). 

[46] G. A. Hagedorn, Commun. Math. Phys. 71, 77 (1980). 

[47] J. Vanderplas and A. Xuereb, “General numeri¬ 
cal solver for the ID time-dependent Schrodinger’s 
equation.” https://github.com/jakevdp/pySchrodinger, 
license: BSD style. 

[48] J. S. Bao, Weizhu and P. Markowich, J. Comput. Phys 
175, 487 (2002). 

[49] H. Metcalf and P. van der Straten, Laser Cooling and 
Trapping, Graduate Texts in Contemporary Physics 
(Springer New York, 1999). 

[50] R. Grimm, M. Weidemuller, and Y. B. Ovchinnikov, 


Adv. Atom. Mol. Opt. Phy. 42, 95 (2000). 

[51] K. Henderson, C. Ryu, C. MacCormick, and M. Boshier, 
New J. Phys. 11, 043030 (2009). 

[52] C. Muldoon, L. Brandt, J. Dong, D. Stuart, E. Brainis, 
M. Himsworth, and A. Kuhn, New J. Phys. 14, 073051 
( 2012 ). 

[53] H. Rauch and S. Werner, Neutron Interferometry: 
Lessons in Experimental Quantum Mechanics, Wave- 
Particle Duality, and Entanglement, Oxford series in 
neutron scattering in condensed matter (Oxford Univer¬ 
sity Press, 2015). 

[54] J. Byrne, Neutrons, Nuclei and Matter: An Exploration 
of the Physics of Slow Neutrons, Dover Books on Physics 
(Dover Publications, 2013). 

[55] M. Utsuro and V. Ignatovich, Handbook of Neutron Op¬ 
tics (Wiley, 2010). 

[56] V. V. Vladimirskii, Sov. Phys. JETP 12, 740 (1961). 

[57] T. Brenner, S. Chesnevskaya, P. Fierlinger, P. Gel- 
tenbort, E. Gutsmiedl, T. Lauer, K. Rezai, J. Rothe, 
T. Zechlau, and R. Zou, Phys. Lett. B 741, 316 (2015). 

[58] V. V. Nesvizhevsky, H. G. Borner, A. K. Petukhov, 
H. Abele, S. Baefiler, F. J. Ruefi, T. Stoferle, A. West- 
phal, A. M. Gagarski, G. A. Petrov, et al. Nature 415, 
297 (2002). 

[59] H. Chen, R. G. Downing, D. F. R. Mildner, W. M. Gib¬ 
son, M. A. Kumakhov, 1. 1. Ponomarev, and M. V. 
Gubarev, Nature 357, 391 (1992). 

[60] R. F. Alvarez-Estrada and M. L. Calvo, in Materials Sci¬ 
ence Forum, Vol. 772 (Trans Tech Publ, 2014) pp. 51-56. 

[61] G. F. Calvo and R. F. Alvarez-Estrada, Nanotech. 15, 
1870 (2004). 

[62] S. Dabagov, Nanosci. Nanotech. Lett. 3, 24 (2011). 

[63] M. L. Calvo, J. Phys. D Appl. Phys. 33, 1666 (2000). 

[64] B. Rohwedder, Phys. Rev. A 65, 043619 (2002). 

[65] S. Longhi, Laser Photonics Rev. 3, 243 (2009). 

[66] D. Dragoman and M. Dragoman, Quantum-Classical 
Analogies, The Frontiers Collection (Springer Berlin Hei¬ 
delberg, 2013). 

[67] P. W. Anderson, Phys. Rev. 109, 1492 (1958). 





